Investigation of the stability of Hopfions in the two-component Ginzburg-Landau 

model 
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We study the stability of Hopfions embedded in the Ginzburg-Landau (GL) model of two oppo- 
sitely charged components. It has been shown by Babaev et al. [Phys. Rev. B 65, 100512 (2002)] 
that this model contains the Faddeev-Skyrme (FS) model, which is known to have topologically sta- 
ble configurations with a given Hopf charge, the so-called Hopfions. Hopfions are typically formed 
from a unit-vector field that points to a fixed direction at spatial infinity and locally forms a knot 
with a soft core. The GL model, however, contains extra fields beyond the unit-vector field of the 
FS model and this can in principle change the fate of topologically non-trivial configurations. We 
investigate the stability of Hopfions in the two-component GL model both analytically (scaling) 
and numerically (first order dissipative dynamics). A number of initial states with different Hopf 
charges are studied; we also consider various different scalar potentials, including a singular one. In 
all the cases studied, we find that the Hopfions tend to shrink into a thin loop that is too close to 
a singular configuration for our numerical methods to investigate. 



PACS numbers: 74.20.De, 47.32. cd 



I. INTRODUCTION 

Topologically stable knots and other vortex-like struc- 
tures have recently received wide interest within con- 
densed matter physics. One reason for this is the con- 
tinuous development of experimental techniques which 
now allow the production of vortices in various types 
of media. These include, e.g., Bose-Einstein conden- 
sate*^, superconductors^^, super-fluid 3 He£i£ and ne- 
matic liquid crystals^. At the same time the increased ca- 
pacity of supercomputers has made it possible to study 
these structures numerically. Thus, in the last decade 
there have been many numerical studies devoted to find- 
ing stable topologically non-trivial configurations in dif- 
ferent physical systems, including, e.g., topological un- 
knots, knots and vortices in the Faddeev-Skyrme (FS) 
model- ' 10 ! 11 ] 12 ! 13 ] 14 ! 15 , vortices in Bose-Einstein conden- 
sate o 16 ' 17 ' 18 , vortex atom lasers in a two-flavor Bose con- 
densate^ and vortices in superconductors^^!, liquid he- 
liu m 22 ' 23 , liquid metallic hydrogen 24 and possibly even in 
neutron stars 2 ^. 

There are various ways in which a vector field in nature 
can support vortices. For example a velocity field can 
form a vortex, e.g., in a hurricane. The position of the 
vortex is in the eye of the hurricane, where the velocity 
is zero. Such a vortex is not stable, because a velocity 
field can continuously change to zero everywhere. 

In some special materials there can exist vector fields, 
associated to spin or other such property, that cannot 
vanish. Then it is possible to have vortices that are 
both non-singular and conserved, the conservation fol- 
lowing from topological reasons. (One example that can 
be produced in a laboratory is the continuous unlocked 



vortex (CUV) in super-fluid 3 He-A£.) The prototype 
model that contains non- vanishing topological structures 
characterized by a Hopf charge was presented by L. Fad- 



dee' 



,9.26 



and recently it has been shown numerically that 



this model does indeed contain stable topological soli- 



tona 10 ! 11 ^ 2 ! 13 . 

In this work we study numerically a system of two 
electro-magnetically coupled, oppositely charged Bose 
condensates described by the Ginzburg-Landau (GL) 
model. Babaev et al. studied this system in Refill and 
argued that it should contain stable knotted vortex soli- 
tons with a non-zero Hopf charge. Here we present our 
results for the GL model, the main conclusion being that 
Hopfions in the GL model seem to be unstable. 

The paper is organized as follows: In Section |TT] we 
formulate the equations, discuss the possible potentials 
and present the initial states that are used in the com- 
putation. In Section IIIII we present the discretization of 
the GL Lagrangian and the numerical method for finding 
stable minimum energy configurations. In Section [PvH we 
describe the results and finally in Section [V] we give some 
concluding remarks on the results obtained. 



II. MODEL 

A. Ginzburg-Landau model 

The model describes two electro-magnetically coupled, 
oppositely charged Bose condensates, as given by the GL 
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Lagrangian density 



Next we define the unit vector field n by 



C 



h 2 

2m i 



2m 2 



(1) 



where we have used SI units. In Eq. © 'J'i and \I/ 2 
are the order parameters for the condensates, A is the 
electromagnetic vector potential, B the magnetic field, 
B=±VxA, and V is a potential. 

Babaev et al£Z- introduced new variables by setting 



* Q = \/2m a px a , 
where the new complex field x is normalized as 



lxi| 2 + IX2| 2 -l, 



and therefore the real field p is given by 



1 /|*i| 2 , l* 2 | 2 



2 \ mi TO2 
In terms of the new fields of ([2]) Eq. © becomes 



(2) 
(3) 

(4) 



V-i£A)x2 



h 2 (Wp) 2 + V(xi,X2,p 2 ) + ^B 



(5) 



The Lagrangian © is invariant under the gauge trans- 
formation 



— i-r- 9(x) .7. 

^ + cd^e(x), 

and the corresponding Noether current is 



(6) 



■Jfc: = |£(**0fc*i-*i0fc*;) 



l*2| 2 



Afc, 



(7) 



mi ni2 
which, using the Ampere's law, must satisfy 

J = V x £? = iV x V x A (8) 
When we use the new variables ([2]) here, we obtain 
J k = 2etip 2 i(x* 1 d k xi - XidkXi ~ X* 2 dkX2 + X2dkX*2) 

c h 

= 4ehp 2 (\j k - f c A k ) , (9) 

which contains a new (non-gauge invariant) current 2 ^ 

3k = i(x*9kXi ~ XidkXi ~ X*2 d kX2 + X2d k X*2)- (10) 
Later on we also use a gauge invariant vector field 

C := jrKj. (11) 



Xi 



X1X2 + X1X2 



n = [Xi X2j o \ I = I i(xiX2 ~ X1X2) > ( 12 ) 
Aj 1 IXi| 2 -|X2| 2 

where a are the Pauli matrices. The inverse transforma- 
tion is 



Xi = n 1 -tn 2 e %a 



X2 = 



_ y/Q--n 3 ) i 



(13) 



V2 



where the phase must be chosen so that Xa are contin- 
uous, see, e.g., (IBTfl) . (This inverse transformation is not 
used in the numerical simulations.) 

In order to show the similarities between the GL and 
FS models, we write the Lagrangian ([5]) in terms of n, p 
and C, after which the n part of the result should be 
similar to the FS model, given by 

£fs = \d k nid k n l + g FS {n ■ d k n x din) 2 . (14) 

After inverting Eqs. © and (fTTj) to obtain A in terms 
of C and j the kinetic part of Eq. © becomes 



^kinetic = fry(|V X l| 2 + I V X 2 | 2 

+ h 2 (v P ) 2 + ^fc 2 . 
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n 



(15) 



Using Eqs. fl2]). © and dTUD one finds that 

<WV = 4(|V X i| 2 + |V X2 | 2 ) - f (16) 

and therefore the first term in Eq. (I15p corresponds to 
the first term in Eq. (fT4|) . 

By direct substitution of Eqs. ([9]) and (fTTj) into -B, we 
also find that 

B = iVxl= A( V x J-iVxC), (17) 

and again using Eqs. (fT2")l . © and (fTU|) . we get 

±e feim n • <9 z r? x <9 m n = -e Mm dij m , (18) 

from which we can see that the B 2 term in Eq. © con- 
tributes to the second term in Eq. (fT4|) . 

Combining the above results we can write the La- 
grangian © in the form 



C =**fd k n l d k n l + h 2 {Vpf + ^fd 2 + V(p, n k ) 



(19) 



which is the form derived by Babaev et al£^. The dynam- 
ical fields are now p, ft and C . If p = constant and C = 0, 
the GL model reduces to the FS model in Eq. (fl~4"|) . Since 
the FS model contains stable topological structures with 
non-trivial Hopf charge, one can hope that the GL model 
also contains similar structures. However, GL contains 
the additional fields p and C in comparison to FS and 
the role of these new fields must be investigated. 
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B. Form of the potential 

A typical and rather general quartic potential used in 
the GL model is 



Vb(*i,* 2 ) = 5Ci|*i| 



\C 2 \*2\ A 



cal^il 2 !^! 2 



+h\V 1 \ 2 + b 2 \V 2 \ 2 + a . (20) 

When ty a are expressed in terms of p and ft and the whole 
system is rescaled so that m a — > 1, we find 



|*i| 2 = p 2 (l + n 3 ), |* 2 | 2 = p 2 (l 
and then the potential (|2U)) reads 



n 3 



(21) 



Vo(p 2 ,n 3 ) = p [n 2 3 {\{ci + c 2 ) - c 3 ) + n 3 (ci - c 2 ) 
+g(ci +c 2 ) + c 3 ] 
+p 2 K(&i-^) + foi+&2]+ao- (22) 

One important aspect in choosing the potential is that 
at infinity the fields will settle to the minimum of the po- 
tential. Furthermore, in order to define the Hopf charge 
it is necessary that the ft field points to the same direction 
far away, otherwise we cannot compactify the 3D-space. 
It would therefore be optimal to have a potential with a 
minimum that would fix the ft field completely, say, to 
n 3 = 1. 

For a particular example assume that c := c± = c 2 = 
c 3 > 0, b := bi = b 2 , then n 3 disappears from Eq. (|22|) 
and the potential minimum is at p 2 = —bj (2c), ft being 
free. For more generic parameter values the extrema are 
obtained for particular values of n 3 and p 2 : 



n 3 



P 2 = 



bi(c 2 + c 3 ) - b 2 {ci + c 3 ) 
bi{c 2 - c 3 ) + 62(01 - c 3 ) ' 
&i(c2 - c 3 ) + 6 2 (ci - c 3 ) 



2(cic 2 - c§) 



(23) 
(24) 



This is a minimum, if cic 2 > c 3 . Note that the above 
values do not necessarily fall within the allowed values 
for p 2 and n 3 , (i.e., p 2 > and |n 3 | < 1), in which case 
the extrema are on the boundaries of the allowed values. 

From physical arguments the following special case is 
relevant 



Vi(*i,* 2 ) =A 



*i| 2 -l) 2 



l* 2 | 2 -l 



(25) 



This breaks 0(3) to 0(2) and corresponds to two inde- 
pendently conserved condensates. It has a minimum at 
n 3 = 0,p 2 = 1. 

Another physically relevant 2 - potential is 



V 2 (*i,* 2 ) 



A ((l^a-i^ + d^l 



1 



-C *i*2 _ *2*1 



an, 



(26) 



which breaks 0(3) completely. This corresponds to 
multi-band superconductors with inter-band Josephson 



effect 29 . In terms of ft this potential is given (using Eq. 

(USD) 

V 2 = A±n 2 + clm[(n 2 - mi)e i2a ] + b . (27) 

The minimu m of this potential is located at p 2 — 
1/2, n 3 = \/l — c 2 /A 2 , while the specific values of n\ 
and n 2 also depend on a, which is related to the phase 
of*. 

From the point of view of Hopf-charge conservation, 
the possibility of p = is problematic, even if it happens 
locally, since then the field n is not defined. However, 
it has been argued^ that quantum effects ensure that 
p 7^ everywhere. In classical field theories, like the 
present one, such expected quantum effects can be in- 
cluded through effective potentials. Since the main pur- 
pose of this effective potential is to guarantee that p 7^ 0, 
its exact form is not so important. One such potential, 
which we will use later, is 



V cfi (*i,* 2 ) = |A(|*!| 2 + |* 



2 2\ 2 
21 - Po) 



27(|*i| 2 + |*2| 2 ) l +ao. 



(28) 



The constants po and a,Q are determined by requiring that 
the minimum of the potential is at p 2 = 1 with a value 
of 0, yielding pi = 2 - 7/A and a = -7 - fa. 



C. Initial states of Hopf invariant Q 

We are interested in the minimum energy configura- 
tions of topologically distinct configurations of the field 
X, or by Hopf-map, ft, related to the complex physical 
fields * Q through Eqs. (J2|) and fl3|). From the point of 
view of ft it is only necessary that limi^i—^ ft = rioo is the 
same in all directions. From this it follows that we can 
compactify R 3 — ► S 3 and then ft becomes a map S 3 — * S 2 
with homotopy group ir 3 (S 2 ) = Z, characterized by the 
Hopf charge. 

Therefore, we have to create a configuration with 
* : M 3 -> C 2 and A : M 3 -> M 3 , such that ft has the 
desired property mentioned above. For * Q this implies 
|*i | 2 + |*2| 2 7^ everywhere and limi£|_ >00 ^ a are (inde- 
pendent) constants. The first condition is also necessary 
for defining the field x an< i if the the second condition is 
also satisfied x becomes a map S 3 — > S 3 . 

The method of constructing a configuration with a de- 
sired Hopf charge has been investigated by Aratyn et 
al£^. They used the fact that for any function <fi ■ S 3 — > 
S 3 combined with the Hopf map h : S 3 — * S 2 




(29) 



the Hopf charge of h o (f> : S 3 — > S 2 equals the degree 
of the map : S 3 — > S 3 (the degree is the n-dimensional 
generalization of the 1-dimensional degree, also known as 
the winding number). 
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For an explicit construction one uses the toroidal co- 
ordinates (r],£,,ip) of M 3 defined by 



sinh(?7) cos(ip) 



A 



X3 



sin(0 



x 2 



sinli(ry) sin((/j) 
A ' 



(30) 



A = cosh(?7) — cos(£). 



Fixing n or Xi does not say anything about the mag- 
nitude of From energy consideration we must choose 
lirri| £|_>oo \I/ Q to be one of the minima of the potential and 
this fixes a preferred value for p. We use this preferred 
value for all x when constructing the initial configura- 
tions. 



In these coordinates the core is at 77 = 00, while the z-axis 
and spatial infinity are at rj = 0. 

Next take any monotonic function 3: [0, 00) — ► [0,1], 
(or [—1,0]) choose p, q E 1 and define the map (j>: S 3 — > 
S 3 by 



0= (g(v)cos(p£),g(7i)sm(p£), 



yl-g(v) 2 cos(g^), \Jl -g(if) 2 wa{q<p)). 



(31) 



If furthermore g is such that g 2 {oo) — <? 2 (0) = ±1, then 
the combined map, ho <p, has the Hopf invariant 

H(ho<f)) = ±pq. (32) 

For details, see Ref^. Using (|3Tj) we now identify 

xi'=<h + i<h = g{v)j*, (33a) 

(33b) 



X2 



At the z-axis X2 looks like a t/?-vortex, and therefore for 
continuity we add the further requirement that 3(0) — 
±1. Similarly, around the core (located at 77 = 00) we 
have a £- vortex, and therefore we demand that 3(00) = 0. 
From Eqs. (T2]) and ((33]) we obtain 



(34) 



Finally, to close the loop, the x of (133]) can be recovered 
from (fTB"]) using ([M]) and choosing a = —qip. 

Note that since 77 = at infinity and g(0) 2 = 1 the 
above construction implies Hoo := limi^i^oo n = (0,0, 1). 

Inverting the toroidal coordinates defined in Eq. ([3171) 
enables us to express Eq. ([3"3"]) in the Cartesian coordi- 
nates. In addition, denoting r 2 = x\- 



2g (v) V 1 - g 2 cos(K 

= ( -23(77) Vi- s'sinCK 
2 5 2 (77) - 1 




ing p = 1, 3(77) = l/cosh(?y) we obtain 



- a; 2 and choos- 



*2(x) 



^/ (r 2_ 1) 2 +te 2 

r 2 + l 
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1 Vv/^f+^f; 



(35a) 
(35b) 



This is the formula used for the initial configurations of 
^ a for our numerical computations. The state indeed has 
the correct Hopf invariant, which, in the case of p = q = 1 
can be seen by finding the preimages of ^ a = 0, which 
correspond to preimages of 77.3 = ±1. These form two 
closed loops, namely r = 1, x 3 — for = and the 
z-axis, which in the compactified space S 3 is actually 
closed. 



III. NUMERICS 

In this section we will describe the discretization of the 
model and the method of the minimization of the La- 
grangian |T]). We will also compare the GL to FS model 
and present some test calculations and the parameters as 
well as coupling constants for the simulations. 

It is important to note that all our simulations are 
done with the fields ^> a only, the fields n, p and C are 
never used in any computation, they are only used in the 
analysis of the results. 

It is perhaps useful to mention once again that the 
change of variables in Eqs. ©-((I]) is not reversible when- 
ever p(x) = 0. There is, however, no guarantee that p 
stays non-zero in numerical simulations of the physical 
fields ^ a and A, even if the minimum of the potential is 
at a non-zero value of p. The situation when p = lo- 
cally can imply breakdown in topology and therefore in 
our numerical simulations we have monitored the changes 
in the global minimum value of p 

Topology can also break if the topological structure 
shrinks smaller than the lattice unit length. In order 
to be aware of this possibility we have monitored the 
global minimum of the dot product of nearest-neighbor 
n-vectors: when the global minimum becomes negative 
the lattice is probably too coarse, or the configuration 
has a genuine singularity. When this has happened in 
the simulations, we have repeated the simulation with 
ever increasing lattice size until the computing resources 
were exhausted. This strongly points to a singular con- 
figuration. 

For simplicity we have used the rescaled (j^ a \ 2 /m a — > 
I 'I'al 2 ) Lagrangian and natural units (c = h = 1) 
throughout our numerical work. With these choices, the 
Lagrangian density used in all our numerical simulations 
becomes 



C 



(V +igl)Vi 



+ §<?/(V x A) 2 + U(* 1 , 2 ). 



(36) 



This can be further scaled by A — > -A, which reveals the 

g ' 

fact that the only relevant parameter is g 2 /gf = 4/ioe 2 . 
The values used in numerical simulations are 3 = 1 
and 3/ G {0.01,1,2,100}; these amount to a particular 
choices of units for po and e. 

In our studies of the FS model, we have always as- 
sumed that limi^i^oo 77,3 = 1, but some of the present 
potentials, e.g., V%, do not have that as a minimum. In 
this case we may assume that limiaji—Kx, n-i — 1, but this 
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can be transformed to 77,3 = 1 by a global rotation in the 
configuration space, which takes n 3 — > n 2 and n 2 — > — n^. 
The same effect can be achieved by using new fields de- 
fined by 

*i = ^(*i + i *5).*a = ^(*i-»*a). ( 37 ) 
changing Vi of (f!i5|) to 

^=A((|*i| 2 + |* a | 2 -2) a 

- (*!* 2 -****) 2 ). (38) 

The different potential terms also contain various pa- 
rameters. For Vo, we have always used 4ci = 4c2 = 2c3 = 
—61 = —62 = clq, which, denoting A = |ci, enables us to 

write Vq — A(|4'i| 2 + |4 V 2 | 2 — 2) 2 . In the numerical sim- 
ulations, we have used potentials Vq, Vi rot and V e s with 
A G {0,1,4, 100,1000}. 



A. Discretization 

The system has been discretized on a cubic rectan- 
gular lattice (indexed as (s,u,v)) with periodic bound- 
ary conditions. Our model can be considered as a two- 
component version of the time-independent, ordinary 
Abelian U(l) Higgs model, which has long since been 
discretized for lattice simulations in quantum field the- 
ory (for example, see Refj 32 ' 33 and the references therein). 
The main point in that context is to discretize the fields 
so that gauge invariance is preserved. From Eq. ([6]) 
we see that if we use the forward discretization of the 
derivatives, the gauge transformation of Af. has to be 
discretized as follows: 

-^i|s,«,i; — > Ai\ s u v + ^ (6 s +i,u,v ~ Qs,u,v)i (39a) 

A-2\s,u,v ~ * ^.2|s,«,u + a (8s,u+l,v — Qs,u,v)i (39b) 
^3\s,u,v > ^3|s,u,u + „(^s,«,« + l — &s,u,1])i (39c) 

where a is the lattice parameter. Thus, Ak should be con- 
sidered as living on the link between two lattice points 
parallel to the coordinate axis k. Combining this dis- 
cretization with the way the $ a transforms in Eq. {6]) it 
follows that the following combination 

-*t\ s+ l, U ,v*l\s+l,u,v ~ ( 4 °) 

where k = ||, is gauge invariant. If we cal- 
culate its continuum limit as a — * by expand- 
ing in a (e.g., *i| a+ i (11lt) = ®i(x + a,y,z) = 
*i(x, y, z) + ad x ^i(x, y, z) + ^a 2 d 2 ^i(x, y, z) + . . .) we 
obtain -a 2 |(9i +ikAi)*i| 2 + 0(a 3 ). For * 2 with At we 
use a similar expression with k — > —k. Finally, for the 



other components of A we use corresponding shifts, as 
illustrated in Eqs. §§§ and (flO]). 

For the discretization of B 2 we use the expression 

(41) 

where, for example, 

-^i2isui? — Ai su +i v A\ suv 

(42) 

— A2, s +l,u,v + -A2,s,u,u- 

The Ffei are gauge invariant under Eq. (|39[) . and the con- 
tinuum limit of Eq. (gT|) with Eq. JUD is -a 2 (V x A) 2 + 
0(a 3 ). 

For the potential we use a 2 V a \ suv ) . The discretized 
Lagrangian is the sum of all the above terms multiplied 
by —a/2, and then the resulting sum has the continuum 
limit of Eq. ((36)) . In practice we use the cubic lattices of 
sizes of 60 3 . . . 960 3 . Of these just the sizes 240 3 , 480 3 and 
720 3 are used in actual simulations, the remaining sizes 
used only to verify the code, results and discretization. 

The minimization of the Lagrangian has been done us- 
ing the steepest descent method and the Fletcher-Reeves 
variant of the conjugate gradient method. The gradi- 
ents needed have been calculated symbolically from the 
discretized Lagrangian. 



B. Initial states 

In order to start the simulations we have to generate 
initial states with specified Hopf invariants. These ini- 
tial configurations have been made using Eq. (|35[) and 
choosing p such that V^oo) = 0. Three different values 
of Hopf invariant were used: 1, 2 and 4. However, this 
leaves A undetermined. The topology of the system is 
unaffected by A, so in principle, we can choose any con- 
figuration for it. We have used several different initial 
configurations for A; for all the results presented here, 
we have used two choices. The first one is defined by 
the condition C = => A = |f J, by Eqs. {10]) and (fTTj) . 
This choice enables us to test the validity of the pro- 
gram by comparing the energies of the initial states with 
those obtained for the pure FS model. The second initial 
configuration has been constructed by solving Ampere's 
law for A with fixed ty a . Note that since we use gra- 
dient methods to minimize the energy, we already have 
encoded Ampere's law in our program: the gradients of 
the Lagrangian with respect to A^ indeed yield Ampere's 
law. 

As a further test of discretization we used the same 
initial configurations in lattices of different sizes and with 
different lattice constants. With lattice sizes of 240 3 and 
above the energies were within 1% of each other. 
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C. Comparison with previous FS studies 

We have also compared the new calculations with those 
presented in Ref. 13 . By using an initial configuration, 
where p = 1, C = and V = 0, the GL model (19]) re- 
duces to the FS model (|14p. In particular, the magnetic 

term becomes E Mb = 12 £p oe * [ e kim {n ■ din x d m n)] 2 , 
while in the FS model we have Et fs = \gFS\ji ' d^n x 

din) 2 . If we set c = K = 1, e = |, [Xq = 1 and 5^5 = ~ 
(the value most commonly used in our earlier work), we 
have Et fs = 2Em b - We created the same initial config- 
uration with our old and new codes and found that the 
energies of the initial states agree to within 1 % with a 
lattice of 120 3 . 

During the calculations we have monitored the topol- 
ogy of the system. As we saw in Sec. Ill B[ the conser- 
vation of topology is not guaranteed whenever p — > 0. 
Therefore we have followed the value of the global min- 
imum of p 2 . In keeping with our earlier work, we have 
also monitored the dot products of ti(x) -n[x+jt) , (where 
p e {(1,0,0), (0,1,0), (0,0,1)}). When the global min- 
imum of this dot product approaches zero or becomes 
even negative it indicates possible breakdown of continu- 
ity and therefore of topology during the simulation. 

We have also calculated the value of the Hopf invari- 
ant directly from the field configuration. Using the basic 
differential geometric result of H [h o ^ = n) = deg x, we 
can even use our computational field variables. Since for 
any map / € C*°°(S 3 ,S 3 ), 

de § / = / ^pfdf A d/ CT A dfP, (43) 

Js 3 

we can compute the value of the Hopf invariant by nu- 
merical integration. The Hopf invariant can also be de- 
termined by visual inspection of the linking numbers of 
the preimages. 



IV. RESULTS AND DISCUSSION 

A necessary condition for any stable static localized 
structure is its stability against scaling. Let us assume, 
for now, that p ^ everywhere. Then the simple scaling 
argument of Derrick^! gives a necessary condition for the 
stability against scaling. Consider the behavior of the 
energy under scaling x — > Ax. For the FS model in 3D 
the volume element of the energy integral behaves like 
A 3 while in the integrand ([14]) the first term scales as 
A -2 and the second one as A~ 4 . In the integral the two 
terms then have opposite scaling behavior and if each 
is positive definite, one cannot scale the energy to zero. 
Furthermore, by the virial theorem the preferred scaling 
is the one where the two terms have equal value. In 
fact one can show 3 ^ that in the FS model the energy is 
bounded from below by the Hopf charge: E > c|Q| 3 / 4 . 
This has also been confirmed numericall y 10 ! 11 ! 12 ! 13 . 



In the present model the scaling stability hinges on the 
behavior of the first and last terms in Eq. (JTSJ) . These 
terms are positive definite and their scaling properties are 
opposite and thus, as long as they are both non-zero, the 
system is stable. The last term, however, equals the term 
B 2 in Eq. |T]) and could, in principle, vanish. Indeed, if 
the C field obtained such a value that 

£kimdiC m = -e k i m n ■ din x d m n, (44) 

then the last term would vanish, rendering the system 
unstable against scaling. One could say that topological 
charge leaks from n to C, so that B — > 

We have observed the above scaling instability in di- 
rect minimizations of all fields. The minimization does 
not exactly follow the route of a uniform scaling, but we 
have observed that the torus-like un-knot shrinks into 
a thin tube, and the value of the integrand in Eq. (|43[) . 
which contributes to the topological charge, concentrates 
tightly in a region near the vortex core. Eventually the 
nearest neighbor values in the corresponding n are anti- 
parallel, which in effect means that the vortex configu- 
ration has shrunk to a size less than the lattice spacing. 
This is accompanied with p approaching zero, and we 
have verified that p approaches zero at the same location 
where the neighboring fVs are anti-parallel. At a place 
where p = the field n can no longer be determined from 
ty, and the topological structure breaks, after which the 
system rapidly goes to a vacuum state. We have ob- 
served this scenario for various different potential forms 
and strengths. 

Two steps of the process of topology breakdown are 
illustrated in Fig. [TJ where we have plotted two distinct 
preimages of n, n 2 — 1 and n 2 = — 1 and the isosurfacc 
ri3 = on which the tubes lie; the colours (gray-scale) of 
the isosurface describe longitudes. Initially (left image) 
the preimages are spaced 90 degrees apart on the isosur- 
face 713 = 0, but during the minimization, on the region 
of icy-plane inside the torus core, the preimages deform 
so that eventually there are no lattice points between 
the preimages of n 2 = 1 and n 2 = — 1. This happens 
even with V e s (in which case p > everywhere) and 
is not accompanied by a shrinking torus; the torus de- 
forms somewhat but the radius of the inner intersection 
of 77.3 = and icy-plane stays approximately constant. 
The two preimages eventually touch each other and when 
they drift apart again (right image), they have become 
reconnected so that the topology is trivial. After this, 
there is nothing to prevent, the energy from dropping to 
zero. 

It should be noted, however, that for any non-trivial 
configuration $ a , Eq. (T4~4"l) violates Ampere's law, Eq. 
©, because by using Eqs. HI]), (T7D and flgj), Eq. (gH) 
implies B = and V x B = J = 0, which contra- 
dicts the topological non-triviality of Thus any non- 
trivial configuration, where Eq. (|44p holds, is unphysical. 
It can then be argued, that although the scaling route 
leads to instability, it does so via non-physical states and 
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FIG. 1: Snapshots of the preimages in the relaxation of the 
(p,q) = (1,2) un-knot, before (left) and just after (right) 
topology breakdown. The tubes correspond to the preim- 
ages of ri2 = ±1. The coloring on the ri3 = isosurface 
corresponds to the longitude of ft. (Colors available on-line.) 



thus does not imply the instability of physically relevant 
states. 

In order to stay within physically relevant states dur- 
ing minimization, we have used the following procedure: 
After every minimization step applied to the 'F fields we 
check whether Ampere's law ([8]) is satisfied. At the begin- 
ning of a simulation the \E' fields are generated from an- 
alytical formulae, such as ([55)1 . which provide the proper 
topological charge. Then Ampere's law is solved for this 
initial state (using conjugate gradients) until is it "suffi- 
ciently accurate". After this initial step we take conju- 
gate gradient iteration steps for all fields until the con- 
vergence criterion for A is no longer satisfied, at which 
point we continue only with A and Ampere's law until it 
is satisfied to the desired accuracy. This is repeated un- 
til the whole iteration has converged. There are various 
methods to determine when the solution is "sufficiently 
accurate" . Since we use a gradient-based method to solve 
Ampere's law, it is natural to determine the accuracy of 
the solution using the gradients. To accomplish this, we 
observe the absolute values of the gradients affecting A 
and note the maximum of these values at each iteration 
and the initial state before any iterations are made. We 
then compare these maxima with the maximum found at 
the initial state; the solution is considered "sufficiently 
accurate" if the ratio of the current maximum to the ini- 
tial maximum is below 0.0001. 

The result of the minimization process is always the 
same: a singular configuration of <F Q , where the region 
contributing to the topological charge has shrunk to a sin- 
gular line. At the same time B field is generated in the 
hole of the torus. This can be seen from Fig. [21 which de- 
scribes states before and just after topology breakdown. 
In particular one can see how the B is formed mostly on 
a smaller ring in the torus hole and how the initial maxi- 
mum of <Fi at origin vanishes while its toroidal minimum 
becomes disk shaped. The topology is broken in exactly 
the same way as in the unconstrained case. 

The exact type of singularity and the process which 
leads to its formation, depends on the potential. For 
Vq and 14ff, the small radius of any initial torus- type 



isosurface of n.3 shrinks without limit, snaring the region 
of topological interest into a singular loop. For Vi rot , the 
process is the same unless the potential is very strong, in 
which case the region where ri2 7^ 0, shrinks to a surface. 

This process is ultimately a scaling instability, albeit 
a different one from the one considered at the beginning 
of this Section. Now, the scaling is not global, but only 
shrink the smaller circle of the torus. The terms involving 
derivatives of ^ Q grow without limit during the shrinking 
process, but only on a loop or a surface, allowing still the 
total energy to decrease. 

In every case the shrinking continues until the discrete- 
ness of the computational lattice eventually breaks the 
topology. We have tried to follow the development of this 
singularity by increasing lattice density. In each case a 
singularity is reached, but in principle the shrinking can 
stop at some still smaller scale determined by the di- 
mensional parameters of the system (such as penetration 
length, strength of the potential, etc.). 



V. CONCLUSIONS 

We have investigated numerically the two-component 
Ginzburg-Landau model using as initial states torus un- 
knots, with non-vanishing Hopf change Q; this brings 
a topological structure into the system. We have used 
different types of potentials depending on both order pa- 
rameters 'Fq of the system. In all cases the initial torus 
tube shrinks into a thin loop and becomes untrackable in 
our computational discrete lattice: the discreteness al- 
lows the topological structure of the system to disappear 
eventually, contrary to the case of the Faddeev-Skyrme 
model. The topological stability in the FS model is due 
to the fact that the kinetic and topological terms are non- 
vanishing and have an opposite behavior in the scaling. 
In the GL model, the term corresponding to the FS topo- 
logical term is the magnetic field term which contributes 
to the unstability of the system. In a direct minimization 
of the fields, the topological charge leaks from n into C 
allowing eventually Derrick-type instability, while in min- 
imizations respecting Ampere's law the process is milder, 
but nevertheless leads to singularity. 

Thus, it seems that the two-component GL model does 
not support stable topological structures having a non- 
trivial conserved Hopf invariant, due to this scaling in- 
stability. It is still possible that in the GL-system the 
natural size of a stable un-knot is much smaller than in 
the FS-model, or that it is only stable for suitably strong 
potentials. It is also possible to stabilize the topological 
structures by adding a suitable term in the Lagrangian. 
One such term has been introduced in Refi^. These are 
questions that we will study further. 
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(c) (d) 

FIG. 2: Snapshots of a Q = 2 system before ((a) and (c)) and just after ((b) and (d)) the topology breakdown. Cross-sections 
of the system have been taken in xy-plane ((a) and (b)) and rrz-plane ((c) and (d)). In all four panels, the upper two plots 
describe the densities |*?i| 2 (left) and | \I/ 2 1 2 (right) on which the phases of the condensates are superimposed as contours spaced 
7r/3 apart; the lower two plots of each panel describe the energy densities ^gj\B\ 2 (left) and -Btotai — \qi\B\ 2 (right). The 
jagged line in the xy-cross-sections corresponds to (f> = 7r; the jaggedness itself is an imaging artifact caused by the discreteness 
of the data used to produce the image. 
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